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■^J" ' Abstract. - We present a new framework for modeling the statistical behavior of both fully developed 

turbulence and short-term dynamics of financial markets based on the nonextensive thermostatistics proposed 
by Tsallis. We also show that intermittency - strong bursts in the energy dissipation or clusters of high 

ON j price volatility - and nonextensivity - anomalous scaling of usually extensive properties like entropy - are 

naturally linked by a single parameter q, from the nonextensive thermostatistics. 



Scaling invariance plays a fundamental role in many natural phenomena and frequently emerges 
■ from some sort of underlying cascade process. A classical example is fully developed homogeneous 

isotropic three-dimensional turbulence, which is characterized by a cascade of kinetic energy from 
large forcing scales to smaller and smaller ones through a hierarchy of eddies. At the end of the 
cascade, the energy dissipates by viscosity, turning into heat. Recently, some authors 1-3 have 
studied the phenomenological relationship among financial market dynamics, scaling behavior and 
hydrodynamic turbulence. Particularly, Ghashghaie et al? conjectured the existence of a temporal 



information cascade similar to the spatial energy cascade found in fully developed turbulence. 

Traditionally, the properties of turbulent flows are studied from the statistics of velocity differences 
v r (x) = v(x) — v(x + r) at different scales r. As with other physical systems that depend on the 
dynamical evolution of a large number of nonlinear ly coupled subsystems, the energy cascade in 
turbulence generates a spatial scaling behavior - power-law behavior with r - of the moments (v™) 
of the probability distribution function (PDF) of v r (the angle brackets () denote the mean value 
of the enclosed quantity). For large values of the Reynolds number, which measures the ratio of 
nonlinear inertial forces to the linear dissipative forces within the fluid, there is a wide separation 
between the scale of energy input (integral scale L) and the viscous dissipation scale (Kolmogorov 
scale 77). Though at large scales (~ L) the PDFs are normally distributed, far from the integral scale 
they are strongly non-Gaussian and display wings fatter than expected for a normal process. This 
is the striking signature of the intermittency phenomenon. After publication of the Kolmogorov 
K62 refined similarity hypotheses 8 , the problem of small scale intermittency became one of the 
central questions on isotropic turbulence. Over the past years several papers 9-16 have discussed 
intermittency and the so-called 'PDF problem'. Similar attempts 1,2 ' 17 have been made to explain 
the same peculiar shape observed in the PDF of price changes z T = z(t) — z(t + r) at small time 
intervals. 
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Based on the scaling properties of multifractals, Tsallis 4-7 has proposed a generalization of Boltzmann- 
Gibbs thermostatistics by introducing a family of generalized nonextensive entropy functionals S q [p] 
with a single parameter q. These functionals reduce to the classical, extensive Boltzmann-Gibbs 
form as q —> 1. Optimizing S q [p] subject to appropriate constraints 6 , we obtain the distribution 

p q (x) = [l-(3(l-q)x 2 ] 1 ^/Z g . (1) 

The appropriate normalization factor, for 1 < q < 3, is given by 
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In the limit of q — > 1, we recover the Gaussian distribution. 

The above distribution, we claim, provides the simplest and most accurate model for handling the 
PDF problem. To show this, we stay in the context of fully developed turbulence (x = v r ). From 
equation (|l|), we can easily obtain the second moment 

and the flatness coefficient (kurtosis) 

^ _ 3(5-3g) 
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We remark that the flatness coefficient, which is directly related to the occurrence of intermittency, 
is solely determined by the parameter q. Also, we note that the positiveness of K r sets an up- 
per bound on the value of q, namely q < 7/5. This bound coincides with the one obtained by 
Boghossian 18 through a ^-generalization of Navier-Stokes equations. Moreover, this limit implies 
that the second moment of distribution (||) will always remain finite, which is empirically expected 
from the phenomena here analyzed. 

At this point, if we assume 2 ' 9,10 a scaling of the moments (v™) of v r as r^ n , the variation with r 
of the PDF of the velocity differences and of its related moments can be completely determined. 
Under the assumptions of Kolmogorov log-normal model 8 , £ n = n/3 — yg/in(n — 3), where \i is 
the intermittency exponent. Then, we immediately obtain the functional forms of the flatness 
coefficient and the parameter q, respectively 

K r = K L {r/L) a (4) 

and 

Q 5-7(r/L)« 

with Kl = 3, the expected value for a Gaussian process, and a = —4/^/9. The correspondent 
expression for (3 can be derived similarly from equation @). However, to account for the well 
known asymmetry of the velocity distributions we may consider (3 = (3 + , for v r > 0, and (3 = f3~ , 
for v r < 0. In this case, both the second and third moments of the modified PDF shall be used to 
determine (3~ and (3 + . Equations (|3|) and (||) remain unchanged, as far as the asymmetry of the 
PDF is small. 



We checked our model with turbulence statistics data taken from reference 2, provided by Chabaud 
et al. 13 . Firstly, we adjusted by least-squares to the experimental data, the parameters of our PDF, 
considering both the symmetric and asymmetric forms of equation (|l]). Then, we compared the 
estimated values of q with the predictions of equation (||), for a intermittency coefficient of /i = 0.25 
(the best estimate currently available is 12 0.25 ± 0.05). The results are displayed in Figs, la and 
2. A good agreement is observed through all spatial scales and for all orders of normalized velocity 
differences. Note that the solid line in Fig. 2 is not simply a fit to the data, since [i is obtained 
independently. We also computed value of q at the Kolmogorov scale and obtained q v ~ 9/7. 

The same approach adopted in turbulence can be straightforwardly applied (with x = z T ) to model 
the statistics of price differences in financial markets, as far as the relevant parameters at the integral 
time scale - time span for which a convergence to a Gaussian process is found - are available. We 
tested our model with price changes data taken from reference 2, provided by Olsen & Associates. 
The results are displayed in Figs, lb and 2. Since we do not have an independent estimate of the 
intermittency coefficient of the information cascade, we also adjusted equation (||) to the data and 
obtained a = —0.17 and tl ~ 2.2 days, the corresponding integral scale of the process. We observe 
that the proposed model reproduces with good accuracy the statistics of price differences over all 
temporal scales. However, we remark that the integral scale value of 2.2 days strongly disagrees 
with other estimates 1 of tl (roughly 1 month). 

Nonextensivity, a matter of speculation in some areas 19 , is an essential feature of Tsallis generalized 
thermostatistics. If we suppose a scenario of a cascade of bifurcations with n levels, and scale the 
generalized entropy S q , averaged over a volume of size V = if and normalized by Boltzmann 
constant, we have at the first level 

S q (2V) = 2S q (V) + (l-q)S 2 q (V) (6) 

and at the top of the cascade 

S q (2 n V) = 2S q (2 n ~ 1 V) + (1 - q)S 2 q (2 n - 1 V) ~ 2 n S q (V) + (1 - q)2 n - 1 S q (V) , (7) 

where the higher order terms in S q have been discharged. Cascade processes are also described 
in terms of fractal or multifractals models 11,20-23 . Within these frameworks, in high Reynolds 
number turbulence, the energy dissipation is not uniformly distributed within the fluid but rather 
concentrated on subsets of non- integer fractal Dp dimension. This picture leads to a scaling 
behavior with dimensionality not equal to that of the embedding space. In this case, if we consider 
the cascade of bifurcations described above, we find 

S q (2 n V) = 2 nDF / 3 S q (V) . (8) 

It follows immediately from equations (^) and (|8|) that 
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Note that the parameter q, through equation (|9j), offers a quantitative picture of the transition 
from small-scale intermittent, nonextensive, fractal behavior to large-scale Gaussian, extensive 
homogeneity. For higher values of n and q = 1, at the top of the cascade, we have Dp = 3. At the 
bottom of the cascade (n = 1), using the value q v ~ 9/7, estimated previously, we get Dp — 2.33, a 
good approximation of the so-called 'magic' value of 2.35, often measured in different experimental 



contexts 24 27 . This result and the divergence of the correlation function of the energy dissipation 
in the fluid 9 
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suggest a more stringent bound on q, namely q < 9/7. If true, this new bound immediately implies 
that, in fully developed turbulence, relation {j]jL) a = 2 is invariant, regardless the Reynolds 
number. 

Though only qualitatively, the above picture may be also applied to the information cascade. 
However, since there is nothing equivalent to viscous damping in the dynamics of speculative 
markets, the information cascade depth is only limited by the minimum time necessary to perform 
a trading transaction. 
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Figure 1 (a) Data points: standardized probability distribution p q (v r ) of velocity differences 
v r (x) = v{x) - v(x + r) for spatial scales r = 0.0073L, 0.0407L, 0.3036L, 0.7150L, with L/rj = 454 
and L and r\ being, respectively, the integral and Kolmogorov scales (see text); data taken from 
ref. 2, provided by Chabaud et a/. 13 ; Solid lines: least-squares fits of modified PDF @); from top 
to bottom: q = 1.26, 1.20, 1.11, 1.08; (3~ = 0.69, 0.66, 0.55, 0.62; /?+ = 0.88, 0.82, 0.76, 0.70 (for 
better visibility the curves have been vertically shifted with respect to each other). 

Figure 1 (b) Data points: standardized probability distribution p q {z T ) of price differences z T = 
z(t) - z(t + r) for temporal scales r = 0.0035r L , 0.0276t l , 0.2210t l , 0.8838t l , with t l = 186265 s 
being the integral scale (see text); data taken from ref. 2, provided by Olsen & Associates; Solid 
lines: least-squares fits of modified PDF (|); from top to bottom: q = 1.35, 1.26, 1.16, 1.11; 
0r = 1.12, 0.83, 0.75, 0.75; 0+ = 0.98, 0.72, 0.61, 0.77. (for better visibility the curves have been 
vertically shifted with respect to each other). 

Figure 2 Dependence of the parameter q on normalized spatial (r/L) and temporal scales (t/tl); 
turbulence data: solid line (model), open squares (data, asymmetric adjust), solid squares (data, 
symmetric adjust); financial data: dotted line (model), open triangles (data, asymmetric adjust), 
solid triangles (data, symmetric adjust). 
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